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Einstein static universes are unstable in generic f(R) models 
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We study Einstein static universes in the context of generic f(R) models. It is shown that 
Einstein static solutions exist for a wide variety of modified gravity models sourced by a barotropic 
perfect fluid with equation of state w — p/p, but these solutions are always unstable to either 
homogeneous or inhomogeneous perturbations. Our general results are in agreement with specific 
models investigated in that past. We also discuss how our techniques can be applied to other 
scenarios in f(R) gravity. 



I. INTRODUCTION 

In 1917 Einstein tried to find a static solution of the 
field equations of general relativity that describes a ho- 
mogeneous and isotropic universe. As static solutions, 
in general, do not exist in this setting, Einstein intro- 
duced the cosmological constant A to make the solution 
static [l|. It was noted soon thereafter by Eddington 
that this solution is unstable with respect to homoge- 
neous and isotropicperturbations. However, subsequent 
work by Harrison Q and Gibbons [1] indicates that the 
issue is more subtle than originally thought. In particular 
Gibbons showed that the Einstein static universe max- 
imizes the entropy for an equation of state with sound 
speed c s > 1/V5- These results have been further inves- 
tigated in H where it was shown that the Einstein static 
universe is neutrally stable against small inhomogeneous 
vector and tensor perturbations, and also neutrally sta- 
ble against inhomogeneous adiabatic scalar density per- 
turbations if the sound speed satisfies c s > l/y/5, and 
unstable otherwise. These results allow for the interest- 
ing scenario where the universe may have started out as 
an Einstein static universe @, 0] , thereby allowing for a 
natural beginning for inflation. 

Because of its interesting stability properties and its 
analytical simplicity, the Einstein cosmos has always 
been of great interest general relativity and its exten- 
sions. These static cosmological models have been con- 
structed in braneworld models [13, EH, ITfl [l7| s tring 
theory (l9j ]. and loop quantum cosmology |20L l2l|. In 
addition, models with non-constant pressure have been 
considered [1, ©, EE EH, 03 • 

The first stability analyses of the Einstein static uni- 
verse in f(R) modified gravity can be found in Refs. [22l . 
[23j . In the latter paper, it was found that stable so- 
lutions do exist which were unstable in general relativ- 
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ity. Subsequent work on the Einstein static universe in 
higher order gravity theories [U, (25| , especially existence 
and stability, had led to slightly differing results. For in- 
stance, in [24j it was noted that that there exists only one 
functional form of f(R) that admits an Einstein static 
universe. On the other hand, in (23| the stability was 
analyzed for a model of a different type. Both findings 
seem to be inconsistent at first glance, however, in this 
paper we are able to reconcile all previous results and 
show their consistency. (Also see Ref. (2(| for a discus- 
sion of the existence of Einstein-static models in more 
general modified gravity scenarios.) 

Although the cosmological constant was soon dis- 
missed after its original introduction in 1917, recent ob- 
servations seem to indicate that the Universe is currently 
undergoing a phase of accelerated expansion, consistent 
with the presence of A (now called dark energy) in the 
Einstein field equations. The idea that dark energy may 
indicate the need for a gravitational theory beyond gen- 
eral relativity has recently inspired a vast amount of re- 
search, in what is know as modified gravity or higher 
order gravity. Such models are not new and have been 
analyzed ever since the field equations have been formu- 
lated in their original form. However, it is only recently 
that these models are investigated in the context of al- 
ternatives to dark energy. f{R) models which have a 
viable cosmology were analyzed, and it was found that 
the models satisfying cosmological and local gravity con- 
straints are practically indistinguishable from the ACDM 
model, at least at the background level [13, HI], for re- 
cent reviews see [H, H(J HH ■ However, such models must 
also be consistent with cosmological structure formation, 
which means the study of perturbation theory in modified 
gravity is necessary [3^, [H, [34[ ■ Therefore, the relatively 
simple Einstein static universe is an ideal test bed for 
perturbation theory to gain insight into the principal dif- 
ferences between general relativity and its modifications. 

We analyze the stability of the Einstein static universe 
against homogeneous and inhomogeneous scalar pertur- 
bations in the context of f(R) gravity. In the following 
section we discuss the existence of the Einstein static 
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universe as a solution of the cosmological field equations 
for two types of models, fine-tuned and non-fined-tuned 
models and make connections with previous work. In £11111 
we analyze the perturbations about the Einstein static 
universe. The perturbations can be characterized by two 
parameters, namely the equation of state w and a pa- 
rameter a which depends on the form of f(R) and the 
matter density. We discuss our results in the final ^IVI 



II. EXISTENCE OF EINSTEIN STATIC 
UNIVERSES IN f(R) GRAVITY 

A. f(R) field equations 

In this paper, we consider modified gravity models gov- 
erned by the action 



S = I d 4 xy / —g 



(i) 



Here, f(R) is an arbitrary function of the Ricci scalar R, 
C m is the Lagrangian density of matter, and k 2 = 8irG. 
The field equations associated with this action are well 
known: 



f'Rab - -^fgab 



[g ab U - V„V 6 )/' = n 2 T ab , (2) 



where □ = V m V m , and we use a prime to denote deriva- 
tives of functions with respect to their arguments; i.e., 
/' = f'(R) = df /dR. Note that the Ricci scalar depends 
on the second derivatives of the metric, so the last term 
on the right-hand side of ^ involves fourth-order deriva- 
tives of g ab . The trace of Eq. (f2J) gives 



f'R - 2/ + 3D/' = n 2 T. 



(3) 



We also assume that the matter content of the model is 
a single perfect fluid 

Tab = (p + p)u a u b +pg ab , P = wp, u a d a =—dr,. (6) 

ao 

In these, w is the constant equation of state parameter 
and ao is the constant radius of the universe. We have 
selected the spacetime coordinates (rj,x,8,ip) to be di- 
mensionless. Putting these assumptions into the f(R) 
field equations yields that the Ricci scalar is equal to a 
constant Rq fixed by the radius 

Ro = 6a 2 , (7) 

and that the values of / and its first derivative at R = Rq 



fa = f(Ro) = 2k 2 P , f Q = f'(Ro) = 



Note that this can be rewritten as 



K 2 pal(l + to) 



(8) 



(9) 



R=R 



In other words, if we view the equation of state parameter 
w as fixed, then any f(R) model that satisfies (HJ) for 
one particular value of the curvature R = Rq admits an 
Einstein static solution [3fJ. For example, if M is some 
mass scale the choice 



f(R) = M 2 exp 



3(1 + w) R 



M 2 



(10) 



does not admit an Einstein static solution unless the cur- 
vature is fine-tuned to R = M 2 . 



This gives a (possibly nonlinear) dynamical equation for 
the Ricci scalar sourced by the trace of the stress energy 
tensor. The existence of such an equation means that 
instead of regarding (fj]) as a fourth order equation for 
g a b, w e can regard ([2]) and as a system of coupled 
second order equations for g ab and R. This observation is 
the basis for the treatment of f(R) gravity as a particular 
type of scalar-tensor theory, see e.g. [35} . 



B. Einstein static solutions 

To find an Einstein static solution of the field equa- 
tions we adopt the metric ansatz 

ds 2 = a 2 {~dr] 2 + -,,d()'d() J ). (4) 

where jij is the metric on the 3-sphere 

j.jde^e 1 = d X 2 + sin 2 X (d0 2 + sin 2 9d(/) 2 ). (5) 



C. GR limit with cosmological constant 

General relativity with a cosmological constant A can 
be recast as an f(R) model with 

f(R) = R-2A. (11) 

In this case, we have an Einstein-static solution for 

Ro = 6(1 + w) k 2 ^ = 2 
A l + 3tu ' A l + 3to' 1 ' 

If the equation of state parameter and cosmological con- 
stant are given, these equations determine Rq and p 
uniquely; i.e., in order to have a static configuration in 
general relativity the radius of the universe must be fine- 
tuned. From this it follows that any fluctuations in the 
universe's radius away from this fined-tuned value will 
result in time-dependent cosmologies, which is what we 
will see in mill 
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D. Models without fine-tuning 



where 



It is possible to construct f(R) models such that one 
can find Einstein static solutions for any choice of the 
radius a , or conversely the curvature i?o- To find these 
models, we regard © as an ordinary differential equation 
that holds for all curvature. The solution is 



/(«) 



(13) 



where the mass scale M is an integration constant. 

Goswami et al. [lij have previously considered this 
case and claimed that this choice of f(R) is the only 
one that admits Einstein static solutions. This state- 
ment is perhaps misleading: (TT3|) is the only f(R) model 
that allows for Einstein static solutions of arbitrary ra- 
dius. In other words, one does not need to fine-tune 
the radius or curvature in these models to have an Ein- 
stein static solution, which is quite different from the 
general relativity case above [33]. We will see in i jHI Gl 
that the spectrum of homogeneous linear perturbations 
of this model admit static solutions, which is not possible 
when f(R) =R-2A. 



X a b — habOfy — g a bh Cd V c Vdfo 

+ (g ab g cd - S a c 5 b d ) { [\V m h cd - W {c h d) m ] V m /£ 

+VV c V d /^ + 2V (c ^V d) /^}. (19) 

We can also linearize the equation of motion ([3]) for R, 
which leads to 

- 3(V a h ac - ±V c fc)V c ^ = K 2 6T\, (20) 

where m is the effective mass of the tp field, which is given 
explicitly by 



m 2 = 



3 



MS 



Jo 



(21) 



Equations (fT5|) and (|2U|) are the main equations governing 
perturbations of generic f(R) models. Notice that they 
are presented as a pair of second order linear equations 
for h a h and (p respectively. Using the basic definition of 
ip (|16|) , it is possible to combine these two equations into 
a single fourth order equation for h ab which reflects the 
fourth order nature of the f(R) field equations. 



III. PERTURBATIONS 



A. Linearized f(R) 



B. Perturbations about constant curvature 
solutions 



We consider fluctuations of the background geometry 
and matter content parameterized by 

i a6 , 5g ab = -h ab , 



Sg a b 



(14) 



where the metric fluctuation is understood to be small 
h a b *C g a b- The variation of the Ricci tensor induced by 
perturbations of the metric is 



SRac 



■\g bd v a v c h M -\g bd v b v d h ac 

+ y bd v b v c h ad + y bd v b v a h cd . (15) 

In all formulae, g a b, Rab, fa, /o, etc- refer to background 
quantities, and all indices are raised and lowered with g a b- 
From this expression, it follows that the perturbation of 
the Ricci scalar is 

ip = 8R= ~h ab R ab + g ab 5R ab 

= -h ab R ab -Uh + V a V b h ab , (16) 

where h = h a a - Note how we have defined the scalar 
quantity cp to be precisely the variation of R. Finally, we 
have the linear variations in / and /': 

(17) 

We now linearize the f(R) field equation (|2|) about the 
background solution 



/o (7W + gabOif - VaVbf) + fo(5Rab 



kfohab + X ab = K 2 STab, (18) 



The Einstein static solution has constant Ricci curva- 
ture, which implies that /q, /q and /q are independent of 
spacetime position. That is, all gradients of these quan- 
tities vanish in equations (TTg|) and (|20[) . These then sim- 
plify to 

fo(Rab<P + gabOlfi - V a V b ^) + 

fo(8R a b - h^gab) - hfohab = ' 5T ab , (22) 



and 



with 



3fU(D - m 2 )<p = k 2 ST\, 



m 



2 _ Jq_ _ R 

3fS 3 ' 



(23) 



(24) 



Notice that these equations also govern perturbations 
about locally de Sitter or anti-de Sitter solutions in arbi- 
trary f(R) models. (Equations of motion for perturba- 
tions of de Sitter backgrounds were first derived in [HI].) 



C. Scalar perturbations 

Working in the longitudinal gauge, we can write scalar 
perturbations of the model as 

hab = 2^u a u b + 2$(g ab + u a u b ), (25a) 
ST a b = S P u a u b + u a D b q + u b D a q + S P V a b , (25b) 
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where V a b and D a are the spatial projection tensor and 
derivative, respectively: 

V a b = S\ + u a u b , D a = V\d h . (26) 

In this gauge, the perturbed metric reads 

ds 2 = al [-(1 - 2^)drf + (1 + 2$)j ij d6 i dd j ] . (27) 

From this, it follows that 'J represents the perturbation 
to the Newtonian potential and $ represents the pertur- 
bation to the spatial curvature. In the matter sector, q is 
related to the perturbation in the fluid's 4-velocity and 
the density and pressure perturbations are rewritten as 



Here we have defined 



Sp = p5, Sp = c 2 s 5p = wpS, 



(28) 



where S is the relative density perturbation and we have 
used that the sound speed c 2 = w for a single perfect 
fluid. 

It is useful to perform a harmonic decomposition of 
these scalar potentials and the perturbation to the Ricci 
scalar ip: 

* = V n (v)Y n (0 i ), * = *n(v)Y n (0 i ), S = S n ( v )Y n (9 i ), 
q = qn(v)Y n {O i ), ip = <Pn{v)Y n {9 1 ) . (29) 

In these expressions, summation over n — 0, 1,2, . . . is 
understood. The harmonic function Y n = Y n (8 z ) satisfies 



(3) AF„ 



n{n + 2), 



(30) 



where ^ 3 ^A is the Laplacian operator on the 3- 
dimensional spatial sections of the model (i.e., associated 
with jij). 

Putting (|23|) into (|2"2")l and (|2"3"|) and performing some 
algebra, we find that ty n , q n and S n are expressible in 
terms of $„ and (p n : 

* =* +£ v a + 

JO « OO 

, _ [(1 + w)f (<p n + fcV) + 4m 2 />»] 
2(l-3H/o/^ ' 

where we have used an overdot to denote derivatives with 
respect to r/; i.e., dip n /dr/ = ip n . Once these are used to 
eliminate \& n , q n and 5 n , we obtain the following equa- 
tions of motion 



= Ax, 



where A is a constant matrix with entries 
6-fc 2 

Ai2 



A 2X = 



A 9 



3 ' 

a(k 2 +3)(1 + w) +3 
9p, 2 (l + a + wa)(l + wY 
2(k 2 - 3)(1 - 3w)(l + a + wa) 



5ak 2 w + 3k 2 w 2 a + 2ak 2 + 6 
3(1 + w)a 



(33a) 
(33b) 

(33c) 
(33d) 
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PL rn 2 



/> 2 \ 3/„# 



K^p 



The solution to the equations of motion are simple and 
given in terms of four constants of integration cj: 

x(jj) = Xi (c 1 e + ^ 1 " + cae-^ 1 ") 

+ x 2 (c3e+^" + c 4 e-^ r ') , (35) 

where Xi and x 2 are eigenvectors of A corresponding to 
eigenvalues — uj\ and — u> 2 , respectively. Explicitly, the 
frequencies are given by 



wi = V -4 + VS, w 2 = V -4 - VB, 



where 



-4 = ^(l + ^)fc 2 -l+ 1 , 
2 (1 + u>)a 



1 



1 



1 — w 



a k z 



+6w - 1 + 



6w 



(1 + iu)a (1 + w) 2 a 2 



(36) 



(37a) 



(37b) 



Obviously, we will have unstable perturbation modes if 
either Im(wi) or Im(w2) are nonzero. This leads to the 
following stability criteria: 

Im(^i) = Im(w 2 ) = & A > and A 2 > B > 0. (38) 

In other words, if the above conditions are met for a given 
choice of (a,w,k), the associated perturbations will be 
stable. It should be stress that the stability of a particu- 
lar Einstein static model is completely determined by w, 
k and the value of f(R) and its first and second deriva- 
tives evaluated on the background solution. That is, it is 
not necessary to know the full functional form of f(R) 
to determine the behavior of perturbations, we just need 
know about the first few terms of the Taylor expansion 
of f(R) about the background curvature. 



D. GR limit 



Notice that since /q = for f(R) = R — 2A, the scalar 
mass m 2 = oo in the GR limit. The corresponding limits 



for the frequencies 



lim lo\ = — 

a^O (1 + U!)a 



arc 



= ±00, 



lim u> 2 — w(k — 3) — 1. 

a— >0 

(39) 

The second of these matches the results of Barrow et al. 
0. The first frequencies are formally infinite, and repre- 
sent an artifact of the reduction of fourth-order to second- 
order gravity. The second set of frequencies implies that 
the model is stable for all w(k 2 — 3) > 1. 
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FIG. 1: Regions of stability in the (a, in) parameter space 
for homogeneous perturbations of Einstein static universes. 
Note that there exist unstable k — modes for all models 
with w > 0. 



E. Homogeneous perturbations 

In the case of homogeneous perturbations we set k = 
in the above expressions. We then find that Einstein- 
static solutions will be stable in two distinct regions of 
parameter space. The first we call the "normal region" : 

normal region = {(a, to) \a < (l + u;) _:t (l-6iu) -1 , 

W > -1,Q > > . 



1 1 

w < 

- 3 2a 



1 - 



16 



1 



(40) 



The term normal comes from the fact that the this re- 
gion has w G (—1,0], so the matter sound speed is sub- 
luminal. The other stability region is called the "phan- 
tom region" , and is given by: 



phantom region : 



{(a, to) |a> (l + w) _1 (l-6w) _1 , 
w<-l,a<0}. (41) 



The stable models in this region have w < — 1 which im- 
plies a superluminal sound speed, hence the term "phan- 
tom". Also, the fact that a < implies that either the 
scalar mass m is imaginary, or the density of matter p is 
negative. Moreover, the effective gravitational coupling 
constant is also negative. In either possibility, models 
in the phantom region are quite strange and we are led 
to conclude that these models should be regarded as un- 
physical. 

In Figure [T] we plot the normal and phantom stability 
regions for homogeneous perturbations. We also show the 



general relativity stability region as inferred from (|39|) . 
The key feature of this plot is that all models with equa- 
tions of state w > are unstable with respect to homo- 
geneous perturbations. 



F. Inhomogeneous perturbations 

We now turn our attention to the behavior of inho- 
mogeneous perturbations k ^ 0. Recall that the spheri- 
cal symmetry of the 3-dimensional spatial sections of the 
model dictates that the value of the wavenumber is dis- 
crete k 2 = n(n + 2). Note that as in general relativity, 
the n = 1 mode corresponds to a gauge degree of free- 
dom related to a global rotation. For the n > 2 modes, 
we find that the model is stable for (a, w) lying within 
two regions: 



right region = 
left region = 

where 

3afc 4 



7i = 



{(a,w) \ w > 7i +72, a > 0} . 
{(a, w) | w > 7i - 72, a < 0} . 

+ 2(3-5a)/c 2 -6(2a + 3) 



6a(fc 4 - 2fc 2 - 6) 



(42a) 
(42b) 

(43a) 



and 



9) fc 4 + 
576a 2 l 



(43b) 



7 2 = (6a)- 2 (fc 4 - 2fc 2 - 6P [9fc 8 a 2 - 

12a (a - 3) fc 6 - 4 (35a 2 + 39a 

24 (4a 2 + 3a - 9) k 2 + 324 - 

In Figure [2] we show the lefthand and righthand stability 
regions for the largest wavelength (non-gauge) inhomoge- 
neous perturbations k 2 — 8. As in Figure [T] we also show 
the relevant stability condition for general relativity. 

From Figure [2] we see that all models with w < are 
unstable to n = 2 perturbations. It is not hard to confirm 
that this it true for all n > 2 perturbations. It is clear 
that the stability regions in Figures [1] and [2] do not over- 
lap. From this we can conclude that it is impossible to 
construct an Einstein static universe in f(R) gravity that 
is stable with respect to both homogeneous and inhomo- 
geneous perturbations, which is the main conclusion of 
this paper. 

G. Non-fine-tuned models 

As an example of the general results we have derived, 
we examine the non-fine-tuned models of £|II Dl Using the 
specific form of f(R) given in (|13[) . we can easily obtain 
a as a function of w for these models: 

l + 3w 

a = -, t-, r- (44) 

(1 + k;)(1-3iu) v ' 

For the homogeneous (fc = 0) perturbations, this leads 
to the following frequencies: 

, 6w 



l + 3u; 



1 ± sgn 



(45) 
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FIG. 2: Regions of stability in the (a, w) parameter space 
for k 2 = 8 (n = 2) perturbations of Einstein static universes. 
Note that all models with w < are unstable with respect to 
these types of perturbations. 



It is clear that for any choice of w, one of these frequen- 
cies is zero and the other is nonzero. This means we 
can always find static and homogeneous perturbations of 
these models. Recall that the unique feature of the non- 
fine-tuned scenario is that there exists Einstein static 
solutions for all radii ao- Hence, these static solutions 
are easy to understand: they represent the deformation 
of one Einstein static solution into another one with a 
different radius. We can also see that the non-zero fre- 
quency will be real only for w € (—1/3, 0]; i.e., the k = 
perturbations will be stable when w S (—1/3,0]. 

When the particular form of a in (|4"4")l is inserted into 
the inhomogeneous stability criteria (|4"2")l , we find that all 
n > 2 perturbations will be stable if 



IV. CONCLUSIONS 

In this paper we showed that one cannot construct 
an f(R) modified gravity theory such that an Einstein 
static universe is stable with respect to both homoge- 
neous and inhomogeneous scalar perturbations. Our re- 
sults and general formulae are valid for any generic form 
of the function f(R). Our approach allows us to ver- 
ify results previously derived based on specific forms of 
f(R), see [23]. Furthermore, we were able to reconcile 
the results of [23[ and [24J, which at first glance seem to 
be contradictory. 

Our results show explicitly that perturbation theory of 
modified gravity theories shows a much richer stability/ 
instability structure than general relativity. This due to 
the fundamentally fourth order nature of the theory, as 
evidenced by the 2x2 matrix equation of motion ([52")) 
for scalar perturbations. This should be compared to 
the equation of motion in general relativity (as in Q , for 
example), which is just a single second-order ODE. In 
case of homogeneous perturbation we found a region in 
the parameter space where the equation of state is in the 
phantom regime w < — 1, see the phantom stable region 
of Fig. [TJ It is unclear how this stability region can be 
interpreted from a physical point of view. 

Finally, we would like to point out how the formulae 
and techniques in this paper can be applied to other sit- 
uations. We note that the equations of ^III Al can be 
applied to any f(R) background solution. The formulae 
of ^IIIBI can be used for any constant curvature solution 
of the modified gravity field equations; i.e., they could 
be used to address the stability of black hole solutions of 
the modified field equations. Our method of determining 
stability criteria entirely in terms of the background val- 
ues of / and its first and second derivatives is applicable 
beyond the Einstein static solution. Indeed, it can be 
an efficient method of analyzing perturbations of given 
constant curvature spacetime in the entire space of f(R) 
models. 



w > 



v/5- 1 
6 ' 



(46) 
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